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Fluctuation and dissipation are by-products of coupling to the ‘environment.’ The Caldeira- 
Leggett model, a successful paradigm of quantum Brownian motion, views the environment as a 
collection of harmonic oscillators linearly coupled to the system. However, symmetry considerations 
may forbid a linear coupling, e.g. for a neutral particle in quantum electrodynamics. We argue that 
absence of linear couplings can lead to a fundamentally different behavior. Specifically, we consider 
a heavy particle quadratically coupled to quantum fluctuations of the bath. In one dimension the 
particle undergoes anomalous diffusion, unfolding as a power-law distribution in space, reminiscent 
of Levy flights. We suggest condensed matter analogs where similar effects may arise. 
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I. INTRODUCTION 

Brownian motion is the prototype of classical stochas¬ 
tic phenomena; its quantum counterpart is more complex 
due to the intricacies of quantum mechanics. Dissipation 
of energy, as well as stochastic fluctuations, are conse¬ 
quences of couplings of the system (e.g. a particle) to a 
reservoir. A paradigm of quantum Brownian motion is 
the Caldeira-Leggett model in which the reservoir is mod¬ 
eled as a collection of harmonic oscillators, an adequately 
realistic prescription at low temperatures [1, 2]. This ap¬ 
proach has been applied to a plethora of systems such 
as the spin-boson problem [3], quantum tunneling [1, 4], 
quantum interference [5], including several platforms in 
solid state physics [6]. The notion of quantum Brown¬ 
ian motion has also been used in other contexts, see e.g. 
Ref. [7]. 

A crucial hypothesis in the Caldeira-Leggett approach 
is the linearity of the interaction between system and 
reservoir [1]. The resulting model is sufficiently simple 
to allow for integrating out the environmental coordi¬ 
nates. Much analytical progress is then possible, mostly 
with functional integral methods [1, 2, 5, 8, 9]. Further 
justifying the above approach, additional nonlinearities 
in the coupling to the environment, which can be treated 
perturbatively [10], can be argued to be irrelevant, in 
the sense that the system is effectively describable by 
linear couplings under certain conditions (see, for exam¬ 
ple, App. C in Ref [4]). Further elaborations include 
an environment composed of interacting fermions [11] or 
spins [12, 13]; the latter [13] demonstrating anomalous 
diffusion of a quantum particle coupled to a spin bath. 
In general, an environment not described by harmonic os¬ 
cillators linearly coupled to the system may give rise to 
qualitatively different behaviors beyond the prototypical 
Caldeira-Leggett model. 


Here, we consider systems in which linear couplings 
to the environment are forbidden by symmetry, leading 
to fundamentally different behaviors due to allowed non¬ 
linear interactions. An important example is a heavy 
(classical) neutral particle, where the absence of charge 
results in quadratic coupling to the electromagnetic fields 
of the environment. We find that such a particle, coupled 
quadratically to the bath, undergoes anomalous diffusion 
at low temperatures, characterized by a power law distri¬ 
bution in space, in sharp contrast to the Gaussian distri¬ 
bution of classical (or linearly coupled quantum) Brow¬ 
nian motion. Confirming the inherently quantum nature 
of this result, we find that usual (Gaussian) Brownian 
motion is restored at finite temperatures for sufficiently 
long times, albeit with a diffusion constant depending 
anomalously on temperature. 

The manuscript is organized as follows. In Sec. II A, 
we introduce the general model under consideration, to¬ 
gether with its equations of motion and action. In 
Sec. IIB, we review previous results for different lin¬ 
ear couplings, before starting the detailed analysis for 
a quadratic coupling in Sec. IIC. While these discussions 
are limited to zero temperature, the case of finite tem¬ 
perature is analyzed in Sec. IID. We conclude in Sec. III. 

II. STOCHASTIC MOTION 
A. Equations of motion and action 

Consider a general framework for a classical particle 
of mass m undergoing dissipative and stochastic motion, 
X{t) denoting its position at time t. The equation of 
motion in an external potential V{X), cast in frequency 
domain, is 

-muj^X{u})+X{dxV}iuj) = F[X,uj], (I) 
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with J^{-} denoting the Fourier transform. The as yet 
unspecified dissipative and fluctuating force F’ is a func- 
tion(al) of ui and X, arising from the microscopic de¬ 
grees of freedom in the environment. Note that a clas¬ 
sical description of the particle is valid only at times 
large compared to intrinsic quantum time scales, such 
as Hfmc^ in the relativistic context considered later {H 
and c denote Planck’s constant and the speed of light, 
respectively). The quantum character of the bath is cap¬ 
tured with appropriate F[X,uj], which can be expanded 
as F[X,u}] = Fo{uj) + x(w)X(a;) -!-•••, keeping the low¬ 
est order terms describing fluctuations and dissipation. 
Here, Fo{oj) is the stochastic force acting on a pinned 
particle, while y(w)X(a;) is a (linearized) deterministic 
force, which may renormalize the intertia m. 

In the absence of the external potential, the equation 
of motion can be organized to first order in X as 

[-muj^ - Xi^)] X{uj) = Fo{uj). ( 2 ) 


Near equilibrium, the fluctuation-dissipation theo¬ 
rem (FDT) relates the fluctuations of the stochas¬ 
tic force Fq{uj) to dissipation by Imx(a;) = 
ta,nh{l3tiu>/2)CFioj), with /3 the inverse temper¬ 
ature, and Cpit) = (1/2) (Fo(t)Fo(O)-I-Fo(0)Fb(0) f^e 
force-force correlator [14]— for useful FDT relations, 
see Appendix A. For linear coupling to an environment 
consisting of harmonic oscillators, Fq is Gaussian 
distributed. Consequently, the distribution function 
of X is also Gaussian, its form entirely captured by 
the two-point correlation function (X(t)X(O)). By 
virtue of the FDT, this implies that the dispersion of 
a particle linearly coupled to the environment can be 
fully constructed from (the imaginary part of) x(w). 
In the absence of linear couplings to the bath, explicit 
computation of the n-point functions (moments) of 
the coordinate X indicates that the behavior can be 
drastically different at low temperatures. 

As a concrete example, consider environment coordi¬ 
nates described by a scalar field <I>(t, a::) in one dimension. 
A convenient way to define the model is via the action 


S' = - dt dx 


4 


+ Si , (3) 


which is simply a continuum description of harmonic os¬ 
cillators. This expression arises in a relativistic quantum 
field theory with light speed c, and may also emerge as an 
effective description in a condensed matter system, with 
c describing the propagation of low-energy modes [15]. 
Furthermore, the particle at position X is locally cou¬ 
pled to the field at (or near) x = X through S/. 


B. Linear coupling 

Let us take the interaction of the particle with the 
environment as 

Si= jdtL^{t,X{t)) , (4) 

where is local and linear in $. In particular, we con¬ 
sider two cases: a) ex <&, and b) (x The 

former describes, for example, the interaction of a (non- 
relativistic) charged particle with the ‘quantum electro¬ 
dynamic’ (QED) field of the bath, while the latter de¬ 
scribes a heavy impurity in a Luttinger liquid where dx^ 
represents the local density [16]. The resulting force is 
Fa. ~ and Fb ~ evaluated at the position of 

the particle. To compute force fluctuations for the par¬ 
ticle at rest, we set A = 0. Since Fo(t) = F[X = 0,t] 
is linear in $, its fluctuations are Gaussian, and can be 
computed directly from Eq. (3). Dissipation can be de¬ 
duced from fluctuations via the FDT as ImXa(a^) a; 
and Imxb(w) for small w as we have shown in Ap¬ 

pendix B. The former produces a typical friction force, 
proportional to velocity /a oc —dtX; the latter giving rise 
to a friction force of Abraham-Lorentz form, /b ~ dfX, 
consistent with Ref. [17]. In both cases, the mean square 
displacement (MSD) of the particle at long times behaves 
as [6, 18, 19] 

lim((A(<)-X(0))2) ~log-, (5) 

T 

where r is a short time scale; see Appendix B for a deriva¬ 
tion. The distribution P{t,X) of the particle’s position 
at time t, starting from X = 0, is a Gaussian with a vari¬ 
ance that expands as logt. While for a general form of 
x(a;) the variance is not always logarithmic [6], a linear 
coupling always leads to a Gaussian distribution function 
in X. 


C. Quadratic coupling 

More generally, the coupling of a system to the bath 
can be expanded as a power series in <I>, generically in¬ 
cluding the first order term. An important exception is 
when a linear term is forbidden on the basis of symmetry. 
A specific example is provided by a neutral (but polariz¬ 
able) particle in the bath corresponding to Eq. (3); this 
may represent the true vacuum with zero-point QED fluc¬ 
tuations, or alternatively the ground state of a condensed 
matter system with an emergent gauge field. The force 
acting on a charged particle is invariant under the simul¬ 
taneous transformation e —>■ —e and $ —>■ — $ where 
e is the charge, and must thus have an expansion of 
the form F ~ “e$” -|- “<I>^” -!-••• (derivatives are not 
shown). For a neutral particle there can be no linear 
term, and the interaction starts at quadratic order in $. 
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Specifically, let us consider a perfectly reflecting parti¬ 
cle such that = 0, i.e., subject to Dirichlet 

boundary conditions, see also Ref. [20]. In one dimen¬ 
sion, this condition cuts the space into two disconnected 
regions which can be treated independently. The energy- 
momentum tensor for the scalar field can be computed 
from the Lagrangian density ^ as 

Tjf = ^ from which the force can be com¬ 

puted as F = Tl{x = 0+) — Tl{x = 0“). With Dirichlet 
boundary conditions, this force takes the form 

F=lidx^Lr-lidx^Rr, (6) 

(again evaluated at the position of the particle) where L 
(i?) denotes the field evaluated at the left (right) of the 
particle. Indeed, the force starts at quadratic order in 
in harmony with our symmetry argument^. 

Employing a second-quantized approach, the field (on 
either side, dropping the L/R subscripts) is expanded as 

$(t,x)= J sin/cx -I-, (7) 

where k is the momentum, and annihilation and creation 
operators satisfy [ 0 ^, 0 ^,] = 2TTS{k — k'). (We have set 
h = c = 1 for convenience.) Note that $ vanishes at 
2 ; = 0 where the particle is located, and the above alge¬ 
bra enforces the canonical commutation relation of the 
field and its conjugate momentum. To compute force 
fluctuations, the following simple diagrammatic represen¬ 
tation is useful: Let us first consider (0|Fo(a;)Fo(—a;)|0) 
where |0) represents vacuum, and Fq(l>j) is the Fourier 
transform of Eq. (6) with X — Q pinned. Fo(—w) is 
a bilinear operator in <I>, and, upon acting on vacuum 
|0), creates two ‘photons’ with momenta k and k' such 
that w = k + k' ] these photons are then annihilated by 
Fo(w) [22], see Fig. 1(a). Each line in the figure is accom¬ 
panied by a factor k which is easily derived by noting the 
derivatives in Eq. (6) and the normalization in Eq. (7). 
One then finds (0|Fo(a;)Eo(—a;)|0) ~ dkk{uj — k) ^ 
for small positive to. An explicit computation for a tra¬ 
jectory X{t) yields the friction force [23-25] (restoring H 
and c), 

Imx(w) = or / = r^X ■ (8) 

OTTC^ DTTC^ 

Computing the second moment of displacement di¬ 
rectly via Eq. (2), or by following the same line of rea¬ 
soning as in case (b) of Sec. IIB, we find a similar result 


for the MSD that depends logarithmically on time, 

((A(t)-X(0))2) = ^log;^. (9) 

A = h/mc is the de Broglie wavelength, and r = h/mc^. 
(We have thus verified that the independently calculated 
Imy and Cp satisfy the FDT.) The force / can be inter¬ 
preted as the friction of vacuum [26] , and is related to the 
dynamic Casimir effect. Since our model is Lorentz in¬ 
variant, / can only depend on derivatives of velocity. The 
two-point fluctuations are thus almost identical to the lin¬ 
ear model (b), and indeed the MSD takes a form similar 
to Eq. (5). Nevertheless, as we shall see, higher-point 
fluctuations drastically change such a correspondence. 

Higher moments of the particle’s position can be com¬ 
puted through Eq. (2); the n-th cumulant is given by 

((A(t) - X(0))")^ = ff[^ - 1) Ri^^) 

x(0|Fo(cc„)...Fo(wi)|0)^. (10) 

The subscript c denotes the connected component of the 
correlation functions^, quantifying deviations from Gaus¬ 
sian behavior [27]; R is the response function obtained as 
the inverse of the bracket in Eq. (2), 


9 / \ • / A ’ 

—rriLU^ — X\^) ~ 'i^sgn[uj) 

with ie inserted to ensure causality. After collecting in¬ 
ertial terms from Rex (redefining m as the total mass), 
higher powers of frequency in real and imaginary parts 
of X (Imx ~ iiof be important compared to the 

leading term at long times. Hence, somewhat surpris¬ 
ingly, xi^) can be dropped from Eq. (11)! 

To go beyond the MSD, we first compute force fluctu¬ 
ations with n = 4 (contributions from L and R cancel 
out for odd n). There are three distinct diagrams con¬ 
tributing to n = 4 as depicted in Fig. 1(b). In contrast 
to the two-point diagram (a), the first two diagrams (i, 
ii) in Fig. 1(b) include vertices with an incoming and 
an outgoing line; the corresponding frequency is now the 
difference of the incoming and outgoing momenta, for 
example, u }2 = k — k' in diagram (i). Thus there is a 
relatively large phase space, k = k', for this frequency to 
vanish. This is not the case for the leftmost (rightmost) 
vertex as both incoming (outgoing) momenta must van¬ 
ish in order to have wi = 0 (a ;4 = 0). The significance of 
this observation is that the integrand in Eq. (10) is quite 
sensitive to the frequency poles in the response functions 


^ One can also formally define an action Sj from which the force 

in Eq. (6) is obtained [21], for example, Sj = a f ^ For example, {{X(t) — XfO))^)^, = {{X(t) — X(0))^} — 3{{X{t) — 

with a 00 . X(0))^)^. 
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Figure 1: Diagrams contributing to the connected n-point 
correlation function in Eq. (10): Each vertex (solid dot) is an 
insertion of Fq with two incoming/outgoing lines representing 
photons, (a) Two-point function; two photons are created 
by the first insertion of Fq and annihilated by the second. 
Energy conservation dictates u = k + k'. (b) Three distinct 
diagrams contribute to four-point functions. Diagrams (i) and 

(ii) have internal vertices with one incoming and one outgoing 
line. The corresponding frequency [for example, the second 
vertex of (i)] generates a pole when k = k'. (c) A typical n- 
point function with all—except the first and the last—vertices 
attached to both an incoming and an outgoing line. 

of Eq. (11). In the vicinity of one such pole, the frequency 
dependence in the first line of Eq. (10) can be cast as 

sm{uJit/2) 1 
uji uji + ie 

Unless there are powers of oji coming from (0|Fo • • • Fo|0)c 
(which is not the case for vertices with one incoming and 
one outgoing line, see Appendix C), this expression is 
divergent near = 0. Therefore, we can use the iden¬ 
tity l/{uji -I- ie) = pA _ iTYS(uJi) to evaluate the inte¬ 
gral. In particular, substituting the delta function in 
the last equation yields a linear dependence on t since 
lim,^;_>o ^ t. Both diagrams (i) and (ii) in 

Fig. 1(b) have two such vertices, and thus contribute to 
the fourth cumulant of the particle’s displacement as 

((A(^)-X(0))4)f’^“^~A2(^;^)^ (12) 

where we have introduced a velocity scale v. In contrast 
to the diagrams (i) and (ii), the contribution of diagram 

(iii) grows as a power of logt. To obtain finite results 
from integrations, we needed to introduce a cutoff fre¬ 
quency, r, beyond which the particle is assumed to be 
completely transparent. The velocity v is related to this 
cutoff by f/c = 2r /m 1; the inequality following from 
the reasonable assumption that T <C to. Equation (12) 
implies a rapid increase with time of the fourth moment, 
compared to the slow logarithmic growth of the MSD. 
In striking contrast to a purely Gaussian process, where 
{{X{t) — X(0))")^ = 0 for n > 2, the fourth cumulant 
grows faster than the variance at long times, rendering 
the process more and more non-Gaussian as time evolves. 

Given the above, we have computed all cumulants of 
the particle’s displacement to the leading order at long 


times. Such computation requires the sum over all dia¬ 
grams with any (even) number of vertices in which each 
vertex, except the first and the last one, is attached to 
one incoming and one outgoing line, see Fig. 1(c). The 
result for the sum of all such diagrams (with A = X/y/Qn) 
is 

2a" 

((A(t)-A(0))")^=—(u^)"-^ (13) 

to leading order in t, and for all (even) n > 2 (see Ap¬ 
pendix C for details). All cumulants (n > 2) grow as a 
power-law in t; the exponent n — 2 coming from the num¬ 
ber of vertices with both incoming and outgoing lines. 
Again, this is in sharp contrast to traditional Brownian 
motion where the cumulants for n > 3 are subdominant 
to the variance. As a consistency check, we remark that 
the moments M„ = ((A(t) — A(0))") satisfy the con¬ 
straint X M„+2m that simply follows from 

positivity of probabilities (Pawula theorem [27]). 

The long-time dominance of higher cumulants in 
Eq. (13) indicates that the probability distribution func¬ 
tion is broad; indeed an infinity of experiments would be 
needed to measure finite moments. Any experimental ob¬ 
servations will be sampled from the distribution function 
P{t, X) for the position of the particle at time t, starting 
from X = 0 at < = 0. Such distribution is consistent with 
the simple form 

= for \<\X\<vt. (14) 

The region jXj ;< A requires a quantum mechanical treat¬ 
ment of the particle; for |A| > rit ^ A, the distribution 
function P is possibly nonzero, but suppressed at long 
times. Remarkably, the distribution function in Eq. (14) 
is a power law, in contrast to the Gaussian distribution 
of the usual Brownian motion (or even quantum disper¬ 
sion with linear coupling). The absence of linear coupling 
to the environment thus qualitatively changes quantum 
Brownian motion. A notable feature of Eq. (14) is the 
emergence of a wavefront that propagates with the speed 
V. The probability distribution is peaked around |A| ~ A, 
and its only time dependence resides in the linear prop¬ 
agation of the wavefront in the tails. Nevertheless, prob¬ 
ability is conserved, f dX P{t,X) « 1, since for ^ A 
the contribution from the tail is negligible. 

Power-law tails of the distribution are reminiscent of 
Levy flights. However, unlike the latter our process is not 
Markovian. In particular, Eq. (14) does not satisfy the 
Chapman-Kolmogorov equation that marks Markovian 
processes [27]. This reflects the long time correlations 
(memory or coherence) of the quantum bath at zero tem¬ 
perature. Introducing a finite temperature indeed makes 
the system Markovian at sufficiently long times. 
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D. Finite temperature 

All the above models acquire a finite mobility at fi¬ 
nite temperature. This may be expected since thermal 
photons colliding with the particle mimic a more classi¬ 
cal Brownian motion. At a temperature T, the loss of 
coherence occurs over a time scale t ~ h/ksT, leading 
to a Markovian process with the standard diffusion law: 
The MSD increases linearly with t, — Ar(0))^) = 

2D(T) t with a hnite, temperature-dependent, diffusion 
constant D(T) that is related to mobility fi{T) via the 
Einstein relation D{T) = kBTfj,{T). The distribution 
becomes more and more Gaussian over time, with cu- 
mulants that are all linear in t, and thus subdominant 
at long times^. Nevertheless, we stress that the anoma¬ 
lous diffusion in the absence of linear coupling [Eq. (14)] 
persists even at finite temperature for times t < h/ksT. 

Here, we consider in some detail the model in Sec¬ 
tion lie at a finite temperature T. It is useful to 
first compute the two-point field correlator. Defining 
<I>'(t) = 0)), we find (in units that ks = h = c = 1) 

poo 

($'(t)$'(0)) ~ / dkk [e-*'=‘(n(fc) + 1) + e*'=‘n(fc)] 

Jo 

^T^csch^intT), (15) 

where n{k) = [exp(A:/T) — 1]”^ is the Bose- 
Einstein factor. In the zero-temperature limit, 
limT->.o (4’'(t)4>'(0)) ^ 1/t^, indicating long-time 

memory, and thus non-Markovianity. The force-force 
correlator at finite temperature follows from Eq. (15) as 

(Fo(t)Fo(0)),-T4csch4(7rtr). (16) 

Equations (15) and (16) show that the correlation func¬ 
tions are exponentially suppressed at long times. Let us 
recall the definition Ci?(f) = (Fo(t)Fo(0)-|-Fo(0)Fo(t))/2. 
We are interested in the long-time asymptotics, and thus 
expand the Fourier transform of Eq. (16) at small fre¬ 
quencies to obtain Cf{uj) ^ + O (w^). The response 

function is related to fluctuations via the FDT. Specifi¬ 
cally, at high temperatures, Imy ^ {ui/T)CF{to), hence, 
we obtain 

lmx(a;)^r"a;, (17) 

leading to a friction force linear in velocity. It then fol¬ 
lows that the mobility scales as 

KT) ~ T-2 , (18) 

giving rise to an anomalous diffusion constant D{T) ^ 
1/T as T —>■ 0. The divergence of the mobility as T —>■ 0 


® This is because, for example, ((X(f) — X(0))^)^ ~ t <SC 

((X(f) — X(0))^)^ ~ at long times. 


is expected since the model is Lorentz invariant in the ab¬ 
sence of a background ‘photon’ bath. Including n-point 
functions beyond n = 2 does not change the above con¬ 
clusions as shown in Appendix D. 


III. DISCUSSION AND OUTLOOK 

Fluctuation and dissipation in a quantum bath are in- 
triguingly different from their classical counterparts. We 
have shown that non-linear coupling to environment co¬ 
ordinates, robustly dictated by symmetry considerations 
in certain settings, can lead to yet another distinction: 
A non-Gaussian power-law distribution (in one dimen¬ 
sion) akin to a Levy flight (albeit non-Markovian). Fur¬ 
ther studies in higher dimensions should be of interest 
from both fundamental and practical perspectives; we 
briefly remark that the argument that led to the power- 
law dependence on time, namely the frequency poles in 
Eq. (10), appears to hold in all dimensions. 

In this work, we have focused on a particle interacting 
with a scalar field in vacuum; however, similar consid¬ 
erations apply to a wide range of problems. First we 
remark that the Lorentz symmetry can be emergent 
(as in Luttinger liquid), while the ‘particle’ may be 
an impurity in a condensed matter system. The key 
ingredient is a general symmetry principle that forbids a 
linear coupling to the environment. As another example, 
consider the transverse-field Ising model at its critical 
point, described by 4)^ theory [15]. For a non-magnetic 
impurity that couples symmetrically to f and ), spins, the 
force should be symmetric under <I> —>■ — $, and start at 
quadratic order. Neglecting the <I>'^ term, such impurity 
undergoes the same power-law diffusion described above, 
assuming it obeys Eq. (1). The <I>'^ interaction, however, 
may change this behavior in subtle ways with different 
universal properties, making for an interesting topic of 
future study. 
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Appendix A: Fluctuation-dissipation theorem 

In this appendix, we provide some useful forms of the 
FDT. The response function R defined by 

(A(f)) - (A(0)) = I Rit- , (Al) 
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is related to the correlation function, C{t) = 
{l/2){X(t)X{0) + X{0)X(t)), via the FDT as [14] 

lmR{u}) = h~^ tanh{/3hu}/2)C{u}). (A2) 

At zero temperature, this equation takes the form 

Imi?(a;) = h~^sgn{uj)C{uj). (A3) 

By computing the inverse Fourier transforms, one finds 
at zero temperature [19] 


^(A(t)-A(O))") 


h r°° 

t' t' Ft t' — t 


where P denotes the Cauchy principle value. In a typi¬ 
cal situation, the response R{t) defined in Eq. (Al) ap¬ 
proaches a nonzero constant value for times large com¬ 
pared to a relaxation time r, i.e., R{t) = fj, for t ^ t, 
which upon the substitution in Eq. (A4) yields the result 
given in Eq. (5). 

We can interchange the role of the dynamical variable 
and the force in Eq. (Al). The equivalents of Eqs. (A2) 
and (A3) are given by 

Imx(w) = tanh(/3fia;/2)(7_F(a;), (A5) 

at finite temperature, and 

Imx(w) = n~^sgn{uj)CF(uj ), (A6) 

at zero temperature. 


Appendix B: Linear models: Dissipation and 
displacement 

In this appendix, we present a derivation of the results 
obtained in Sec. IIB. Consider the quantum field defined 
by the action in Eq. (3). The field in free space can be 
written in a second-quantized basis as (S = c = 1) 

$(t, x) = J ^ , 

(Bl) 

where ujk = |A:|. 

For a particle weakly coupled to the field $ via Eq. (4) 
with oc $, the force is Fo(t) ^ i9a;$(t,0). Some alge¬ 
bra yields 


Fo{t) ~ i 



dk 

27r 



(cfc — a-k) 


e*'=‘(a 


t 

k 


«-fe) 

(B2) 


In particular, for w > 0, Fq{uj) = Fi[(—w) ~ i^/uj (a^j — 
a-uj) where the subscript on a denotes momentum. One 
can then easily see that C'i?(a;) ~ ui; we have used the 
fact that (Ola/jcJ,, |0) = 2TrS{k — k'), while all other two- 
point correlators vanish. Equation (A6) then dictates 
Imx(a;) ^ w. The response function in real time is given 
by R{t) ^ 1 for t ^ t, with some r > 0. 

On the other hand, if the particle is coupled to the 
field via (x the force becomes Fo(t) ~ i9^<i)(t, 0). 
Then, 


Fo{t) 


dk 

3 27r -^2 L 


^—ikt 


(ofc + a-k) + e* (al + a_/.) 

(B3) 


and Fo(uj) = Fq{—oj) ~ (a^ + a-^j) for w > 0. 

Hence, Cpiai) which in turn implies, via Eq. (A6), 

that Imx(a;) ~ The response function (neglecting 
runaway solutions) becomes R{t) ~ t for t 3> t with 
some r > 0. The linear term gives a ballistic propaga¬ 
tion corresponding to the second time derivative in the 
equation of motion. 


Both response functions computed above yield a loga¬ 
rithmic growth of the MSD when plugged in Eq. (A4). 


Appendix C: Nonlinear model: All moments of 
displacement 


In this appendix, we compute higher moments of the 
distribution function defined in Sec. IIC. We consider the 
region a: > 0 with the corresponding force 
(up to an unimportant negative sign); the contribution 
of the region a; < 0 gives rise to an overall factor of 2. 
It is useful to explicitly write 4)'(t) = 9x4’(t, 0) starting 
from Eq. (7), 


$'(t) = ^ . (Cl) 

An insertion of the field produces two incoming/outgoing 
lines, each of which picks up a factor of with k the 
corresponding momentum. Each line is shared by two 
vertices thereby contributing a factor of 2k. Combined 
with 1/2 in the definition of the force, we can assign 
each line a factor of k. However, each vertex gives an 
additional factor of 2 due to the choice of contracting the 
lines attached to it with the incoming/outgoing lines. We 
also write the exponential factors in Eq. (10) in the form 
of sine functions, thus getting an additional factor of 2i 
per vertex. 

These considerations yield 
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((X(t)-X(0)n = 2x 


2^{2iy 


n 

.i=l 


duJi sin(uJit/2) 
2 -k ujf + iesgn{uJi) 


Tn ({w*}) 27r5(wi H-h UJn) , 


(C2) 


where the overall factor of 2 takes into account the contri¬ 
bution from both left and right half of the space, the delta 
function makes the energy conservation explicit, and 
should be computed by summing all u-point force corre¬ 
lation functions. Following the discussion in Sec. IIC, 
we take Wi « 0 for 2 < i < n — I, and a;„ « — wi; this 
allows us to treat wi ^ 1/t as a fast variable and the rest 
of frequencies as slow variables uJi < \/t. The first and 
the last sine functions in Eq. (C2) can be then combined, 
by averaging over the fast variable Wi, as 

sin(a;it/2) sin(a;jit/2) —)• —icos(wt/2), 

where we used conservation of energy, and defined uj = 
W 2 + ■ • • + . Furthermore, the function can be ap¬ 

proximated as J-n{uJiiO, ■ ■ ■ ,0, —uJi) = Fn(iOi). To com¬ 
pute the diagrams in the form of Fig. 1(c), notice that 
there are only two momenta running within each dia¬ 
gram: k and uji — k; this is due to the assumption that 
Wi « 0 for all vertices except the first and the last ones. 
We should sum over all such diagrams with each (except 
the hrst and the last) vertex connecting to lines with the 
same, but either of the two, momenta. It is then straight¬ 
forward to show that 

JL , .n+l 

Moj,) = j -k)[k + {u, - fc)]"-' = 

(C3) 

where the bracket has simply organized the sum over all 
possible diagrams (an expansion of the bracket generates 
all such diagrams). The integration over wi now gives 


^ , F"-2 

io 2it 2ATT-^{n-2y 


(C4) 


where the denominator is due to the response functions 
(due to the first and the last vertices), and F(<C m) is 
a frequency cutoff on the idealized Dirichlet boundary 
conditions, that is, the point particle is transparent for 
frequencies w > F. Finally, the integral over all other 
frequencies together with cos{uit/2) is computed as 


trigonometric function cos{uit/2), and used 



duj sin(a;t/2) cos{ojt/2) 
2 tt + iesgn{uj) 
dw sin^(a;t/2) t 

27r 4 


it 

4’ 


(C6) 


Putting all these together, we obtain Eq. (13). 


Appendix D: Finite temperature 

In this appendix, we extend the results obtained in 
Sec. IID at hnite temperature to higher-n-point functions 
using an approximate, simplified approach. At long times 
{t > H/ksT), we approximately have (ks = h = c = 1) 

Ofe « 4 « yn(k) « ^, (Dl) 

and, Eq. (15) becomes 

/ dk 

^cos(fct) =rd(t). (D2) 

At long times, the dynamics is dominated by the fric¬ 

tional term computed in Eq. (17), 

T^X = F{t). (D3) 

It follows from Eq. (D2) that 

{F{t)F{t'))c T^d{0)S{t - t') - T^S{t - t'), (D4) 

where d(0) is defined by an integral over frequencies 
less than T, hence it is proportional to T. Combining 
Eqs. (D3) and (D4), we find that 

((A(t)-A(O))^) ~ ^ |Y‘dt'dt"(F(t')F(0)c~ 

(D5) 

consistent with the fact that D{T) ~ 1/T. In order to 
compute the fourth moment, for example, first note that 



duJi sin(wit/2) 
27r ujf + iesgn{ujy 

in^n-2 

2”-i 


COS [{uj2 H-1- a;„_i)t/2] 

(C5) 


In computing this integral, we have expanded the 


{F{t)F{t')F{t")F{t'”)), 

- T^5{t - - t'"). 

It then follows that 


(D6) 


(D7) 
















Similarly, we find that all connected moments of X grow 
linearly with time at long times. 
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